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ABSTRACT 


Bie Method of @raruizal images has been successfully 
eaeplted CoO electrostatic peoolceisulnvOolVIne seonduccors on 
a dielectric substrate. This same method is investigated 
for its adaptation to dynamic problems. A Green's Function 
is derived and applied to the problem of a wire dipole 
antenna on a dielectric substrate. The input admittance 
of the antenna is computed by the Method of Moments. 
Experimentally measured values of input admittance are 
aed Wine cm CiCOrcelCalevalues and the crror as 


discussed. 
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LT. INTRODUCTION 


if a monopole antenna is mounted on a perfectly 
eonducting ground plane, it will reflect totally above the 
plane and will produce an image below the plane. This 
imareswiliebave, an image @eurrent in the same direction 
as the monopole in order to satisfy boundary conditions 
SUuMLOCmerouUMmeGnplane. in essence, because of its image, 
the monopole will act as a dipole in a free-space. Hence, 
if the ground plane is removed, the image may be replaced 
by another physically real monopole thereby creating a dipole. 
In other words, conditions above the reflective plane remain 
unchanged if the plane is removed and a real object replaces 


the image. 


A. THE USE OF IMAGE THEORY 

With the monopole, he thickness of the ground plane 
eo SSuUned Lemneeineible, ~ This resules an only one image. 
However, If stmeserouna plane 1S a drelectric with finite 
meeKnessSmeveua, reihectzon mo longer occurs. tinstead, 
mages Pare produced from partial reflection at the surface, 
reflection within the dielectric and transmission through 
the dielectric. These images are used in the "method of 
partial images" [1]. When the dielectric region is removed, 
these images may be considered real charges and are summed 
as an Warinice Series. However, a finite number of them 


gisye™cood results. 





B. A DYNAMIC SOLUTION 

BHecenclNGCmetmparuzaleamares has been successiully 
eppiiied te static problems. The*™ purpose of this research 
is women Vveseleave ats application in solving dynamic 
problems. A searen for a Green's Function iS conducted. 
Once found, it is utilized in an impedance equation. 
Themecguatton weneravres an impedance matrix so that the 
"method of moments" can be used to solve for the current 
distribution or the input impedance of an antenna mounted 


on a dielectric substrate. 





ieee Ne lUne On THE PROBLEM 


Pincee tie Cle@reroOsvaute proolemmuias been Solved, a 
method for solving the dynamic Beep WS sought using 
erina pee neorvemrme dynamic situation issehosen to be a 
half-wave dipole antenna on a dielectric substrate of 
arbitrary thickness and relative permittivity (e,,). A 
wire constitutes the antenna whose radius and length are 


variable. 


A. ON THE IMAGE COEFFICIENT METHOD 
Consider one small segment -of the wire as an element 
of charge. Assume the charge to be a short distance from 


the dielectric. The charge element will have many lines - 


Charge element 


f 






Dielectric, . 


S 


Om Slux in every radial direction. However, in considering 
a single ite i SONnCmen  limpenier rates tne dlelectric, 
Tio ONnCMOtmuaicmret loevedmat the interface. The amount 
reflected is proportional to the reflection coefficient (K). 
Similarly, the part penetrating the dielectric is propor- 
tional to (1-K). 

Naturally, many reflections and transmissions occur 
Sie wvOnvdmomeali tl lemmlixebine. “In turn, it produces many 


image points. The diagram serves to illustrate: 
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(Figure eee) 


Each of@tme three regions, 1,11,1lI1, may have different 
Gielectric constants. However, in this case, regions I and. 
IIT have we (free space). Region II has an arbitrary E 
iecoolvines iermethe Green's Function, each region has a 


different image representation. 


(a). Image representation for left region - I. 
(b). Image representation for center region - II. 


(c). Image representation for right region - III. 
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(i-)e KA-K)q KCI KD 


K?(-K)g -k(I-k)q 


7 (1-k)q K7(1-K)a K*(-K) g 


ENE ge (c) 
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(Figure 2.) 
B. SOLVING FOR THE GREEN'S FUNCTION 
ppecitically, in che probiem at hand) the wire vactually 
lies on the dielectric at the interface of regions II and 
III. Region II is the desired region with which to work 
because the accumulation of charge is assumed to be greater 
Slembine Glielectric Side Of The wire than on the free space 


Side. The image representation for this region appears as 


I 


=H 
-p(l-p)q — ~pli- Pq (i-?) q Pli- P)q Pli- Pq 


where the reflection coefficient is now p= nee in order 
l+er 
to avoid confusion with the propagation constant (k). 


Adding terms to simplify: 


Tce, Vo) WAR 





(b). For We. e(t-p)-e(i-e) = (A P) (p*- #9 
2 Pay = = eiGerns 


’ 


Co). For Was Pie) Pe Cr-P) = (1-0) (P*- 0”) 
= P(1-P)(P-1) -P°(1-e)? 


al 


(api Fer Us, : Pw i ee)= p* (\-0) 
= e'(1-p)(P-1) 


Gee em 
3 2. 
-P (1-9) 
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~ In terms of a summation, the Green's Function for the 


dielectric region is expressed (2-1) 
Lo- 


Wee py -a-er/ Pe Yi 


where Y, = exe(-ily (n-a)* +(y-y')? +(z-2iD)* ) 
And (A-n')? +(y-y')? +(2-24D)? 


V,= ex p(-jk/ (4-2)? + (y-y')4 
An] (t-n)° + (y-y')? 


O = (~€, oe 2% 
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Yet, with a Single wire, variation occurs in eyallity meets 
X-direction. Y 


She 


18, 





The dielectric Green's Function then becomes 


a : i-l) 
Vez (y-e)% - (1-¢) ) is Wi 


where Y. exp (dk f(x-n')@+ (240)* ) 
dar] (t~a')* + (2id)* 


exp (- ak} (a-2')* ) 
Are} (n-ai)* 


il 


Y, 


C. PHASE CONSIDERATION 

in formulating the Green's Function for the static 
problem, the phase of each image point need not be consid- 
ered. However, in the dynamic problem, image point phase 
must be examined. With the charge element a distance "a" 
from the dielectric in Figure 1, each image point has a 
phase term associated with it. This means that Va eourlG 
not be expressed as a Simple summation. 

On the other hand, phase terms are eliminated if the 
ania element ls adjacent Go the dielectric. since the 
Teco memocLUA Ly wadjoined wo the divelectri¢, every image 
Pes0b appearsceasmtnough it were immersed “in the dielectric 


region. For instance, consider image point "q," which 





aly 





CoM Uwe oOlnet mace path Cla this Cl comes from flux line A, 
which is reflected into B and then ©. The image path length 
MOI as equal tothe addition of dielectric path lengths A 
and B through simple geometry. in other words, the image 
path "C1" appears to be in the dielectric medium due to A 
and B. Then image point "ay" also appears to originate in 
the dielectric for the same reason. This is true for every 
image point. Hence, all image points appear to be in the 
same medium with equal phase. The Green's Function, then, 


PolLaLNs 2cs Simple form. 


D. THE USE OF ¥, 

Min Creed is ubbye (ial eye (Y) ef this=secuvien was derived 
in a manner similar to Silvester's solution [1] using the 
MechodsOtm partiabeimages. | Silvester, of course, solved a 
static problem. Yet, it seems logical to assume that his 
method should be applicable to a dynamic problem through 
upe reasons Stated. 

| lic ne CiedmetemOnenivs I1Ssused to findwasnumerical 

solution to the input impedance of the antenna. This method 
utilizes the. Green's Function in solving for the impedance 


Mee 1 dbsce 
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Git. THE EMPRDANCE MATRIX 


fiieoolvi ew Orme ecutercny GitStrabuLion @r the Preed— 
Peril eMpecanees Ole an anvenna, the mavrix AyEROaws in the 
method of moments iS a good approximation. The key is to 
properly VWoadgpne impedanceymatrix. Once this isi accom- 
Plishned, assumed values of voltage fill the voltage eto 


eames che problem is solved. 
2 -1 
Litt) Ge [hvala 


The impedance matrix is obtained through the impedance 
equation. The equation is derived by Harrington [2] ana 
the details are given in Appendix A. The Green's Function 


(Y5) moon UNempMevIOUSsseCClLONn 1S CMploeved in Che Seouacion. 


Z (mn) 7 jopAds inte Wing) ty ta Ries 


AWE 


- Win ih) ; Wa m) + Yee] 


A. CALCULATION OF THE GREEN'S FUNCTION 

In order to express the Green's Function in matrix form, 
Vg must be integrated along the wire. Holding the source 
point (n) fixed while the observation point (m) ranges over 


every sepmenvy generates the ne column Gite the Mea trios Liem 


us: 





MOVingeeche SOurcee point to tite (a) Se Pe venerated. 
allowing the observation point (m) to roam along the wire 
generates the next column, etc. Finally, the (m x n) 
matrix is filled, where "m" is chosen numerically equal 
be eV! 


In terms of equations: 


V (nm) = AU Wa df 
A Xn 
A&n 
| -jkRom n an ikRam 
Ve mn) = (i-p)e@ dy’ - \ (\-py" 0 rant 
Abn A. ag Rom A3Yn A ax Rim 


hat 


where 
= NY 
Rea (hyo ) + 
men 
Re nt [ (zi)? 2)0 
L_™ yA _ _ a) _ 
E> | Chm ) - (a. oa! = (H,. A.) 
mén 


Ri = [ld ma! yee (2e0)? = / (a -a,)%4 (250) 
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There is a different Green's Function for each case. 
The derivations are detailed in Appendix B while the results 
are shown here. 


For case 1. (m=n) 


4 nh es) Sey lyre 
Von m)= (i-p) Onto) ik = (1-9)" a e _han 
Rac | Atm, = R) Aa Z (24D) 


heeecaccmil i. — (mz) 


ik (Amin) 
Vn m) =(\-e) e ; 
4 ax Cain eah) 
a: 7 AK f Caen a)? + (QED 
-(\-?) Aa e° 





A = ae ee 
i iz (Aenean rE 


B. CALCULATION OF THE IMPEDANCE MATRIX 
mie IMpe@amecavallesmare noWwselOund Using ene calculated 


Green's Functions in the impedance equation. 


Z. (m,n)? 4p AQn Von mn) ste Karas 
Hwe 


- Waa) - Vorym) + Vc, my) 
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Where At, = Adn- iN Since the length of each segment 
Hemchiosen to bpescoual, 

Three different situations arise in calculating elements 
Ceeenie amMpedancermavriux ss (ne first is ¢ahe "se1f—impedance" 
term. Whenever the observation point and source point 
Pome ldewmamoe !l—imMpedance term is poner ed) ine se  Wwerns 
occupy the diagonal of the matrix and are all numerically 
equal. 

The second situation is the "mutual-adjacent" term. 
nescmoccurmsawhen the observation point is adjacent to the 
source point. Due to the fact that the distance between 
"m" and "n" is always the length of one segment CAVES) - 
every mutual-adjacent term is equal. 

The third situation happens when the observation and 
source points are one or more segments apart and are called 
"mutual impedance" terms. They complete the matrix and 


very ln Magsnavude. 
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ie elk impedance 
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4+ gt 


] |! a ) (20) 
Ztnyr)= jopAte | aa af 6 (0) we 
Ro 


ri OID) 








+2 gg 





ds hot Ak (23D) 
jk (0) - as 
jo€, | 2a] An 2] Aa (2;D) 
BLAS r p-\) -aK . L . 
-2 es = = Fey f se jk (A Qn) + (2 p) 
AWE, |4ar Ady A at = J CAkn)* + (22D)* 


ei 





2. Mutual-Adjacent Iimpecanee 


- . 
ral ; 


mi yn 


— =p ma -= 


= 
3! 3¢ 


L~t iki (A n\4 iL 
Z (m, Als = jah tn (\- (i-P) enh ee: ) e4 Ona (2i DP 
Aw ASn (2 (Adn)*4 (2: Tbs apy 


jk An f tie - 4k J(A%n)* + (23d)? 
S 


+2 I(1-e) - (i-e) 


ee ein, ed 


{WE An Aldn Ac TE [(Adny + a ))* 


(2i-1) eat 2AM) + (240) 


~\ i (-e) _. ctl 
AWe,| A 2A, a [An + Rid D)? 
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yee el 20) 
= |r| Date) ju] cae Tet 
4We ZW At 21 ter a ais 








3. Mutual Impedance 


i. i Pe 


“ nm om m Clee | ate [ae 


{M= | mn -~j 


cL AQn (21-1) al (LAL2.D) 
Z. (rn, la jou Ale (i-e) e) saat \1-e) De e' ic 
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fe DNs I Afwes (aud) nh2iD) 


ee (Ret) ik (LAR 4 ONG 
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ge, A. ax Yay CAT Ps Qo sae weed 
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{WE Lx (LP A Qn) 7 


KPa (rit) Gk /(LPAGn? + (2iD) 
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‘an PA “TALPA Qn + (2A D)* + Os One 


_ A _ os Neae ne) 
-| {(i-p) e ee 
ha 
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C. PROGRAMMING 

The impedance equations of the previous section are 
programmed for computer usage. Of course, similarities 
exist throughout the equations making the programming 
easier. 

Once the proper parameters and constants are selected, 
the impedance matrix is filled and then inverted. The 
values of voltage in the voltage matrix are chosen to be 
zero everywhere except at the feedpoint, where the source 


is one volt. 
-1 
[Z] {v] = (T] 


The resulting matrix [I] will be the current distribution 
on the wire. The middle value of [I], the feedpoint value, 
is the feedpoint admittance since the source was one volt. 
Tac computer puoeram is written to solve for various values 


Pere COpOINnGemormeanpUG. addmittances. 
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LVR LA VALU A LION 


After programming the equations, meaningful data is 
sought for analysis. Several parameters are available: 
Chem@radius, the Mlength, the dielectric constant, the 
thickness, and the frequency. In order to evaluate the 
validity of the equations, a cemErgiaon is made against 
Harrington's solutions. Hence pie mOskeLOO Ui emeonst any 
iene gion Li ise chosen to be El. The length is equal to 
half of a free space wavelength at f = 3 GHz. The length- 
to-diameter ratio, L/2a, is set at 74.2, while the frequency 
moe varied. 

The results are plotted and compare favorably to 
Harrington's graphs. Figures 3 and 4 show conductance and 
Susceptance, respectively, versus the length-to-wavelength 
ratios. It was found that at least a forty-one by forty-one 


matrix was necessary to achieve accuracy. 


A. EXPERIMENTAL PROCEDURE 

Experimentally, a thin copper strip was attached to a 
one-eighth inch thick dielectric substrate with Ee 1k The 
Gielectric was mounted perpendicular onto the ground plane 
so that the thin strip would be fed as a quarter-wave 


monopole vice a half-wave dipole. The antenna was driven 


luarrington, Pore lod Compuuan lon oy Moment. Methods , 
p. 72, The MacMillan Company, 1968. 


eal 
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Input Admittance of the Dipole 
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through 50 ohm miniature coax and input impedance was 
measured with a network analyzer. 

The experimental values of input conductance and sus- 
ceptance are plotted along with values calculated using the 
impedance equations. Obviously, Figures 5 and 6 illustrate 
a significant difference between experimental and calculated 
results. There is an amplitude variation and a difference 


in resonant frequencies. 
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V. CONCLUSIONS 


trOnemacs Pree=spaceseurvesewhich mateh Harrington's 
curves, it seems that the Green's Function and the 
impedance equations are valid. of counse, Jn this instance, 
e,=l so that the reflection coefficient (9) is zero and 
every term containing (Pp) is eliminated. The part remaining 
is simply a free-space Green's Fucntion. Hence, the solution 
is valid for a half-wave dipole antenna in free-space. 


DieeorilmavelveeneeSare is not true when ea us 


A. EXPERIMENTAL DIFFERENCES 

Aside from the normal attenuation losses of experimental 
measurements, it seems the only other inaccuracy would come 
from the adhesive used to attach the copper strip to the 
Gierecrric.  =tt!igs would Influence the effective dielectric 


Pons t anu | © ) and, perhaps, the amplitude of the admit- 


BFF 
tance. Therefore, the emphasis is switched to the 


calculation procedures. 


B. PROBLEM IN THE GREEN'S FUNCTION 

The ER Lewellen PacmGreeniselunclten was chosen 
alrows! SOMe@error to exist The accumulation of charge 
was assumed to be greater on the dielectric side of the 
wire rather than the free-space side. Because of this, 
Region II was selected in accordance with the method of 
partial images to solve for the necessary Green's Function. 


In the process, the propagation constant is dependent on 
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the relative permittivity off the dielectric. 
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~- 2 


Af€r 





Seer 
Ad 


LOWeVCi en COmen ne DiySmeammesLAandpOlnG. aS tLhewlimiv oF 
the dielectric thickness approaches zero, the Green's 
Me vlna mioulcde: CAUCCHTO =a shree-space soluvion. Yet, W, 
from section TTP will never reduce to a free-space function 
due to the dependence of the propagation constant on E, 
This suggests the use of both a free-space propagation 
constant and a dielectric propagation constant. But there 
does not seem to be any apparent method in combining the two 
Propagation constants*with a Green's Function for Region If. 

Further investigation led to a Green's Function for 
Region Ill whieh incorporated both propagation constants. 

It reduces to the free-space function as zero dielectric 
LALeCkKness 1s "approached. |) However, une soluvion leads to 
negative conductance. 

Finally, Binide the method of partial images can produce 
valid solutions for static and free-space dynamic problems, 
it would soon that with some additional changes this same 
method should solve all dynamic problems. Unfortunately, 
the final answer lies beyond the research efforts of this 


paper. The method of partial images, however, warrants the 


oWeXeNel jone ie bar elevene Tai <eksy eal ekehwalionale 
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APPENDIX A 


DERIVATION OF THE IMPEDANCE EQUATION 


ihe mmpedancesequauszonsis developed in 4a manner similar 


to Harrington's derivation. Formulating the problem: 


E 7 -jwoA- VO (A- 1) 


iad aE | 
A ? uD Js e’ 4s (A 2) 








Aa PR 
~4kR 
= 1% (ar e* dS (A - 3) 
‘ AaR 
OF a an i/. ie (A-4) 
4 
nN xX E* = =n "se on S (a- 5) 





non A (Figure A-1) 


ani 





The wire is broken into many segments subject to the 
BOmLOw Live Cond One . 


(1) Current flow is assumed to be along the wire axis. 


(2) To the axial component of E at the surface, apply 
boundary condition (A-5). 


With these assumptions, equations (A-1) through (A-4) 


become 


ra me 


~Ey =-ywAg- 2& oS (A-b) 
df 
_ -4kR | 
A=p\Tde! df (A-7) 
mis eR 
~jkR 
@=1\el(fh ec af (A-8) 
< ATR 
AtS 
OG = - he dT (A-9) 


a 
= 


J 


= 


i 


\ 


where "{" is the "length variable" along the axis of the 


wire. 
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Figure A-1 illustrates the divistion of the wire into 


"n" segments. 


a + 
Swarts aul, so Centered abn, and terminates av n. 


Recording to thegnotation, any me segment 


The 


distance between any fi and n is denoted by Af, the 


length of the segment. 


Summing integrals over many small 


segments of Ad, approximates an integral over the entire 


Spe yee a 


approximated 


~ Em ~ jw A gion 7 ® (i) 3 D (m) 


KH 
4+ 
rz; 


3+ 
r 


 ( 
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€ 


— | 


Equations (A-6) through (A-9) are then further 


(A-t0) 
Nee 


_ -4KR 
LimAln |e. dh Cat) 
A Ly Am R 
BS 
-jkR 
+ 
Lh CoC) er i (A~\2) 
A&R 
n ARS 
Lays EG (A- 13) 


AX 


ey 





Similarly 


mln) & seer Gel) (A- 14) 
jo A Qs 


substituting equation (A-13) into (A-12), and equation 


(A-14) into an equation similar to (A-12) for ® (ad, 


+ gk 
d (m) = — \ Lat - Lm) e d Q (A-\5) 
jure Ati “hak 
Axn 
-4kR 
@ (rm) = -_ tL Jn) - Lind e a (A-16) 
fe LQ Aka Aw R 
Axa 


halo Om CMe evel Chin VOmCXDreSSethe anverrals parts of the 


above equations in the form of Green's Functions. 


. 


i KR 
Von m) = NO e Af 


= Sees 


. AMS, A wR 
Ay 


Where Ki. ViCwGrsvanice Biron vine midpoint of Uhe nm, 


sepment Go whe anvegration point on the n— 
Serement. 


Using this definition, equations (A-11), (A-15) and 


(A-16) become 


A (rs) = b> it (nn) AY, Vin m) (A-17) 
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® (mh) = ate) [ten eh - Teh, a (A-\8) 


AWE 
Bm) = ne oe Tony Va | (A-\7) 
AGIE 


Now substituting equations (A-17), (A-18) and (A-19) into 


equation (A-10) and noting that A ¢ (1) z Alm) « AYen 


fm 


- Eg (m) = - Gwe im Lind A fae A Pm Ycnm) 
A Men 


flat > LT (nad Y Ga chy ~ Lon) Won mm) 
Abn {we} Z 
aA). Ton Ya, m) ~ L(n- Va mn) | | (A-20) 
Paar 


Two further approximations are made 
Tine Voi my = Ton Vom (A-21) 
and 


Ton) Yaa) = LW Ya a) (a-22) 


With these two substituted, equation (A-20) is written as 
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where 


Z (rn) = jop Ada AL Vem) + 1 | Yee 
Yb 6 8) (A-24) 


Equation (A-24) is the impedance equation which applies 
for the case of (m=n) self impedance, as well as for (m#n) 
mutual impedance. 

Since "m' roams along all "n" segments of the wire, 


there is a set of "n" linear equations. In matrix form: 
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If CV] and [2] are known, the current distribution may be 


© UT ee 
[I] = (2)77 Cv] 


as well as the feedpoint impedance. 
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APPENDIX B 


DERIVATION OF THE GREEN'S FUNCTION 


There is a distinct Green's Function for each case: 
when the observation point coincides with the source point 
(m=n); and when the observation point and source do not 
coincide (mn). 
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Using a Maclaurin series expansion on the first term: 
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Since the firsey two verms ofthis Maclaurin expansion give 


reasonably good results, 
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b. Integration is now performed on the second term of 


the Green's Function. 
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